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Abstract 

The low pt direct photon production from a variety of the hadronic sources is studied within the 
microscopic HSD transport approach for p-\-C,p + Pb and Pb + Pb collisions at 160 A GeV. The direct 
photon emission from elementary hadronic scatterings as well as meson-meson bremsstrahlung are in- 
corporated. The influence of in-medium effects such as a collisional broadening of the vector-meson 
spectral functions on the photon emission rate is found to be hardly observable in the final spectra 
which are dominated by bremsstrahlung type processes. The uncertainties in the subtraction of the 
'background' from the photon decay of hadronic resonances inside the hot and dense fireball is inves- 
tigated, additionally. Our findings are relevant for the interpretation and extraction of experimental 
data on direct photon production at low pt- 

PACS numbers: 25.75-q, 25.75.Cj, 24.10.Lx 
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I. INTRODUCTION 



The properties of hadronic matter under extreme densities and temperatures and their phase 
transition to the deconfined and strongly interacting Quark-Gluon-Plasma (QGP) are the cen- 
tral topics of modern high-energy physics. In order to understand the dynamics and relevant 
scales of this transition laboratory experiments under controlled conditions are presently per- 
formed with ultra-relativistic nucleus-nucleus collisions. The electromagnetic radiation (real 
and virtual photons, i.e. dilepton pairs) is the unique probe to study the heavy-ion collisions 
since the photons do not suffer from the final state interactions with the surrounding matter 
and, thus, provide a clear signal of the various environments of their creations (cf. the pioneer- 
ing; works . Moreover, the photons are emitted from all stages of the collisions - from 
the QGP to the hadronic phase, so they may provide information about the different degrees 
of freedom - from partons to hadrons jj, [|], m, U.\- 

On the quark-gluon level the main processes for the direct photon production are Compton 
scattering of quarks and gluons qg — > q^y and annihilation of quark and anti-quarks qq — > gj as 
leading processes, while the next leading order (NLO) process is dominated by bremsstrahlung 
Q9 Q91- The photons radiated from the different sources have a different transverse mo- 
mentum pt that allows to separate the contributions: The photons with high p? are primarily 
produced by the initial hard NN collisions and called as 'prompt' (or 'hard') photons. Direct 
photons from a thermalised QGP dominate at lower transwerse momentum (1 < pt < 3 GeV/c) 
and denoted as 'thermal photons'. On the hadronic level the main source of direct photons is 
meson-meson rescattering: up — > 7T7, titt — > fry, ttK — > K*j, Kp —>■ Kj, where the first two 
channels are the most important. These photons from the hadronic rescattering dominate at 
even lower px- Additionally there are photons coming from the hadron decays (such as tt°, rj 
etc.), which give a huge photonic 'background' and make the experimental measurement of 
direct photons very complicated. 

Most of the theoretical predictions for direct photons are based on the local thermalization 
assumption and evaluate the photon production rates from the equilibrated quark-gluon plasma 
or hadronic matter which then are convoluted with the space-time evolution of the system. In 
order to account for the non-equilibrium dynamics with its full complexity one needs microscopic 
transport models. However, most of the transport models - commonly used for the description 
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of high energy heavy-ion collisions - are based on the hadron-string picture (e.g. HSD 1 and 
UrQMD 2 ) and do not include the phase transition from partonic to hadronic matter in a 
consistent way, nor solve the hadronization problem. We mention that there are attempts to 
develop such transport models, e.g. the multi-phase transport model (AMPT) [8] which includes 
pQCD-like partonic scattering, or combined models - hydrodynamics for the QGP + UrQMD 
for the hadronic stage or Parton-Hadron-String-Dynamics (PHSD) model which is based on 
a dynamical quasiparticle model matched to reproduce lattice QCD results in thermodynamic 



equilibrium 



io| 



We recall that already in 1995-1997 there were attempts do describe the first experimental 



photon data by WA98 Collaboration on S+Au at 200 A GeV 



ll| in the transport models 



UrQMD Q, AMPT Q and HSD lX The transport calculations included the photon 



production by the hadronic decays jj, Q as well as direct photon production from the 
meson-meson scattering ITa [13fl. Also the influence of in- medium effects such as dropping 

nn 

vector meson masses on the photon yield have been addressed in Refs. |13l . I14J ] . Unfortunately, 
the first WA98 data provided only an upper bound for the direct photon production yield and 
did not allow to draw solid conclusions - all transport results were just below the upper limit 
given by WA98. In 2000 the WA98 Collaboration provided new data for Pb+Pb at 160 A GeV 



15| | . This stimulated a new wave of interest for direct photons from the theoretical side (cf. 
the review [la]) and references therein). 

The aim of our present study is to investigate the photon production at SPS energies from 



hadronic sources using the extended version of the HSD transport model [171 . Il8| which includes 
the off-shell dynamics of the vector mesons with dynamical spectral functions. It allows to study 
the influence of in-medium effects - such as a collisional broadening - on the photon emission 
rate. We will compare our transport calculations to the experimental data for Pb+Pb at 160 A 
GeV as well as to the new preliminary data on photon production in p+C and p+Pb collisions at 
160 A GeV [l9j. We stress here again that the HSD model doesn't include the phase transition 
from QGP to hadronic matter; we thus concentrate on the low photon production which is 
dominated by hadronic sourses. First 'pilot' HSD results have been reported in Ref. j^ . 
Our paper is organized as follows: In Section II we describe the treatment of various channels 



1 Hadron-String-Dynamics transport approach 

2 Ultra-relativistic-Quantum-Molccular-Dynamics 
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for the photon production in HSD. Section III contains a comparison of the HSD results with 
the WA98 data for Pb+Pb at 160 A GeV, while Section IV contains our results for photon 
production in p+C and p+Pb collisions at 160 A GeV. A summary closes this work in Section 
V. 



II. PHOTON PRODUCTION IN HSD 



Our analysis is carried out within the HSD transport model [14|, l2lj, [22] - based on covariant 



self energies for the baryons 23] - that has been used for the description of pA and AA collisions 
from SIS to RHIC energies. We recall that in the HSD approach nucleons, A's, N*(1440), 
N*(1535), A, £ and S* hyperons, S's, S*'s and fi's as well as their antiparticles are included on 
the baryonic side whereas the 0~ and 1~ octet states are incorporated in the mesonic sector. 
Inelastic baryon-baryon (and meson-baryon) collisions with energies above y/s th ~ 2.6 GeV 



(and T/sth — 2.3 GeV) are described by the FRITIOF string model 24j whereas low energy 
hadron-hadron collisions are modeled in line with experimental cross sections. Note that the 
HSD transport approach includes the off-shell dynamics of vector mesons explicitly; for the 
details we address the reader to Ref. 

We consider the following hadronic sources of photon production: 

I. The photon production by a mesonic decays (Tr°,T],r)',u),4>,ai) where the mesons are 
produced firstly in baryon-baryon (BB), meson-baryon (mB) or meson-meson (mm) collisions. 
The photon production from the mesonic decays represents a 'background' for the search of the 
direct photons, however, this background is very large relative to the expected direct photon 
signal. Moreover, there are a severe experimental difficulties in subtracting the photons from 
hadronic decays. For the present study we consider the contributions from the photon decay 
of the following mesons: 

7T°^7 + 7, (1) 

r) -> 7 + 7, 
v' ~> P + 7, 

U — > 7T° + 7, 

ai — >• 7r + 7. 



The decay probability is calculated according to the corresponding branching ratios taken from 



PDG [25|]. The broad resonances in the initial/final state are treated in line with to their 
in-medium spectral functions which will be illustrated in the next subsections. 

II. Additionally to the resonance/meson decay channels the photons can be produced directly 
in elementary collisions of particles. For the present study we consider the direct photon 
production by the scattering processes 

tttt -> pry, (2) 

Tip -> 7T7, (3) 

accounting for all possible charge combinations. We note, that we discard 7-production in 
rescattering processes with strange mesons (such as tt+K* — > K+j, tt-\-K — > K*+j, K+K* — ► 
7r + 7 etc.), since at SPS energies the strange meson density is subdominant. However, as 
indicated in Ref. 26j such processes may play a role at higher energies. In the next subsections 
we will define the contributions of the processes tttt — > pj, up — ► ir~f explicitly. 

III. The photon production from meson-meson bremsstrahlung might play an important role, 
too, as has been pointed out firstly by Haglin and later on by Liu and Rapp [28]. These 
authors have calculated the thermal photon emission rate from the proceses m + m — > m + m+'y 
and found that the bremsstrahlung radiation is one of the dominant channels for low energy 
photon production. 

In the next subsections we define explicitly the treatment of the dominant hadronic process 
for direct photon production. 

A. Direct photon production in tttt — > /ry reactions 

Since at SPS energies the pion density is relatively high, the tttt annihilation channel has to 
be accounted for in the direct photon production. With increasing energies (and pion densities) 
the importance of this channel grows accordingly. 

The basic form of the cross section for the tttt — > /ry reaction has been adopted from Kapusta 
et al. in Ref. [29| where the high-energy photon production has been evaluated in the hot 
hadron gas model. However, in Ref. [29|] the final p meson has been considered on-shell, i.e. 
at the pole mass M = 0.77 GeV. Since in HSD the vector mesons are treated with the full 



off-shell spectral function, which depends on density and momentum 



171 ]. some modification of 



the formulae from Ref. 29[ has to be done in order to account for the broad mass distribution 
of the final p meson. The simplest way is to 'fold' the cross section from Ref. 29[] with the 
mass and density dependent spectral function of the p-mesons. 

We note, that for the present study we use the same form of spectral functions for the vector 
m e 80I1 s a 8 mod e,ed in Re, Q. In p articula , the special fusion of a vector meson V with 



mass M at nucleon density p^ is taken in the Breit-Wigner form: 

2 M 2 T* v (M,p N ) 



A V (M }PN ) = C 



1 • 



(4) 



vr (M 2 -M * 2 (p^)) 2 + (Mr^(M, P7V )) 2 ' 
with the normalization condition for any p^: 

j A v (M,p N )dM = l, (5) 

where Mu m = 2 GeV is chosen as an upper limit for the numerical integration. The lower 
limit of the vacuum p spectral function corresponds to the 2ir decay M min = 2m n , whereas 
for the in-medium collisional broadening case M min = 2m e — > with m e denoting the electron 
mass. Mq is the pole mass of the vector meson spectral function which is Mq(pn = 0) = M 
in vacuum, however, shifted in the medium for the dropping mass scenario according to the 
Hatsuda and Lee 30] or Brown/Rho scaling 3jJ. Futhermore, the vector meson width has been 
implemented as: 



r* v (M, \p\,p N ) = T v (M)+T coll (M, \p\,p N ). 



(6) 



Here T v (M) is the total width of the vector mesons (V = p, to) in the vacuum. For the p meson 
we use 



Tp(M) 



r 



AM) 
(M 2 



r 



MA 2 /<? x 



1 _ 4m 2U/2 



F(M) 



(7) 



go 



(Mq — Ami) 1 / 2 



In 



2 2 
Tq is the vacuum width of the p-meson at the pole mass Mq; F(M) is a form factor taken 



from Ref. 



32j as 



F(M) 



2A 2 + M 2 X 2 
2A 2 + M 2 



(8) 



with a cut-off parameter A = 3.1 GeV. This form factor was introduced in Ref. 32| in order to 
describe the e + e~ experimental data. For the uo meson a constant total vacuum width is used: 
r w = r^(M ), since the to is a narrow resonance in vacuum. 
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The collisional width in (E]) is approximated as 

T co u(M, \p\,p N ) = 7 p N < v a l y N >^ a coU — . (9) 

Po 

Here v — \p\/E\ p, E are the velocity, 3-momentum and energy of the vector meson in the 
rest frame of the nucleon current and 7 2 = 1/(1 — v 2 ). Furthermore, <7y^ is the meson-nucleon 
total cross section. 



As discussed in Ref. Il7| - in order to simplify the actual calculations - the coefficient a co u 
has been extracted in the HSD transport calculations from the vector-meson collision rate in 
A + A reactions as a function of the density p^. The numerical results for r co ^(pjv) then have 
been divided by Pn/Po to fix the coefficient a co u in (jHJ). We obtain a co u ~ 150 MeV for the p 
and a co u ~ 70 MeV for uj mesons which are consistent with the experimental analysis in Ref. 



33] . In this way the average effects of collisional broadening are incorporated and allow for an 
explicit representation of the vector- meson spectral functions versus the nuclear density, p^- 

We mention that also finite temperature effects lead to a sizable broadening of the vector 
mesons spectral functions (also at baryon chemical potential //# = 0). This is essentially due to 
scattering with mesons which may contribute to the total width by 70-80 MeV at a temperature 



of ~ 170 MeV according to the early work by Haglin [34|. But for the present study we use a 
simplified modeling of the collisional broadening width which discard an explicit consideration 
of such temperature effects. However, the temperature effects are partly accounted here due 
to explicit meson-meson interactions which also lead to changes in the vector meson mass 
distribution. Since we find that the spectral broadening of the mesons is practically not visible 
in the final photon spectra, an explicit consideration of temperature effects is beyond the scope 
of this study. 

In order to explore the observable consequences of vector meson mass shifts at finite nuclear 



density - as indicated by the CBELSA-TAPS data 35] for the uj meson - the in-medium 



vector meson pole masses are modeled (optionally) according to the Hatsuda and Lee |30| or 



Brown/ Rho scaling [3l| as 



Mo(Pn) = n ^ M ° , r , (10) 
(1 + ap N /p ) 



where pN is the nuclear density at the resonance decay position r, po = 0.16 fm 3 is the normal 



nuclear density and a ~ 0.16 for the p and a ~ 0.12 for the u meson 33|. The parametrization 



( HOP may be employed also at much higher collision energies (e.g. FAIR and SPS) and one does 
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FIG. 1: (Color online) The 7-production cross section a(s, p) for the 7r _ +7r° — > p~+j, 7r + +7r° — > p + +j 
(upper part - (a)) and 7r + + 7r~ — » p° + 7 (lower part - (b)) reaction for p^ = (solid lines) and in 
the case of the 'collisional broadening' scenario (dashed lines) and the 'dropping mass + collisional 
broadening' scenario (short dashed lines) for nuclear density pn = 3po- The dotted lines correspond 



to the results from Kapusta et al. 



29l | for the fixed p-meson mass M = 0.77 GeV, while all other lines 



show the yield integrated over the p-meson mass M. 



not have to introduce a cut-off density in order to avoid negative pole masses. Note that ffTDD 
is uniquely fixed by the 'customary' expression Mq(pn) w M (l — apN / po) in the low density 
regime. 

Thus, we model the photon-p- meson production cross section in 7T7T reactions in the following 
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way: The total cross section a nn ^ n (s, Pat) is 

<Tw*^rr(s,p N )= j dM — . (11) 

The mass differential cross section is approximated by 

f A(M, p N )dM 

m = ^ n (s, M) ■ A(M, p N ) ■ — , (12) 

/ A(M,p N )dM 

where A(M, pn) denotes the meson spectral function (jlj) for given total width T v ([6]): 



the p meson. In Eq. 



291 ] where the p-meson 



M max = y/s is the maximal kinematically allowed invariant mass o 
f fT2|) a^_+ (s, M) is the vacuum cross section from Kapusta et al. 
mass is considered as a free variable (i.e. not fixed to 0.77 GeV as in |29_j). Thus, formula (fl2"l) 
can be used to model the vector meson production in tttc reactions in the vacuum and in the 
medium, too. 

Fig.[T]shows the 7-production cross section a(s, p) for the 7r~ + 7r — ► p~+7, 7r + +7r° — > p + +7 
(upper part) and 7r + +7r~ — ► p°+7 (lower part) reaction for p^ = (solid lines) and in case of the 
'collisional broadening' scenario (dashed lines) and the 'dropping mass + collisional broadening' 
scenario (short dashed lines) for nuclear density p^ = 3po- The dotted lines correspond to the 
original cross section from Kapusta et al. [29;] for the fixed p-meson mass M = 0.77 GeV, while 
all other lines show the cross section integrated over the mass M. Note that the cross section 



from Kapusta et al. [29( for the fixed p-meson mass diverges at threshold. However, the folding 
over the spectral function of the p meson leads to an endothermic behavior of the cross section 
(cf. the solid line) instead of an exothermic (dotted line). 

As seen from Fig. [1] the cross section in the 'collisional broadening' scenario is basically 
smeared out close to threshold. Only when incorporating additionally a dropping mass the 
thresholds are shifted down in energy such that the production cross sections become enhanced 
in the subthreshold regime with increasing nuclear density p^. 

B. Direct photon production in np — ► irj reactions 

In Fig. [2]we show the production cross section a(s) for the 7r~+p° — > 7r~+7, 7r + +p° — > 7r + +7 
(upper part), n + + p~ — » 7r° + 7, 7r~ + p + — » 7r° + 7 (middle part) and 7r° + p~ — » 7r _ + 7, 
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FIG. 2: (Color online) The production cross section a(s) for the 7r~ + p° — > 7r~ + 7, 7r + + p — ► 7r + + 7 
(upper part - (a)), tt + + p~ — > 7r° + 7, 7r~ + p + — > 7r° + 7 (middle part - (b)) and 7r° + p~ — > 7r~ + 7, 
7T° + p + — > 7r + + 7 (lower part - (c)) reaction. The dashed lines corresponds to the p-meson mass 
M = 0.3 GeV, the solid lines stand for M = 0.77 GeV while the short dashed lines correspond to 
M = 1.0 GeV. 

7T° + p + — > 7r + + 7 (lower part) reactions. Here again we use the vacuum cross sections from 
Kapusta et al. [29] and consider the p-meson mass as a free variable since the production of the 
vector mesons in HSD is realized with respect to the spectral function A(M, p^). The dashed 
lines corresponds to the p-meson mass M = 0.3 GeV, the solid lines stand for M = 0.77 GeV, 
while the short dashed lines correspond to M = 1.0 GeV. 

As follows from Fig. [2] the it + p — » tt + 7 cross section depends very strongly on the mass 
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of the initial p meson - for low M the threshold is shifted to low y/s. Thus, one can expect 
an enhancement of 7 production in case of this channel for the in-medium scenarios due to the 
enhanced population of the low mass p mesons. 



C. Photon production by the decay a± — > ttj 

The photons can be also emitted from the decay a% — > rcy. In spite that the branching 
ratio of this process is not well known experimentally and expected to be very small (we use 
Br{a\ — > 7T7) ~ 1.5 • 10~ 3 [36]), this contribution is also considered in our investigation. The 
study of the a% dynamics is an interesting problem itself since it is related to chiral symmetry 
restoration at high densities and temperatures as pointed out in Refs. [37], |38|, |39|, |40 |. 

The production of di mesons in HSD stems from BB and mB collisions via string excitation 
and decay at high energies and by 7r + p <-> a% reactions at low energies. The mass of the a\ 
meson is distributed according to the Breit-Wigner spectral function - Eq. (jlj) - with a mass 
(p) and density (p^r) dependent width r* . The total ai width T ai can be written as 



r ai(^> \p\,Pn) = T ai (/i,p N ) +T ai . coU (ii, \p\. 



Pn), (13) 



where the collisional width is approximated by T ai:co u(p, \p\,Pn) — oc ai :coii Pn/po with a ai 
150 MeV. 



coll — 



Since the a± meson decays dominantly to tc + p [39( one has to take into account the mass 
distribution of the final p meson A(M,p^) when calculating the width of the a\ meson. Thus 
formula ((7j) for the a\ width T ai (p,p N ) has to be modified accordingly: 



f dM q 3 (p,M) A(M,p N ) 
r ai (p,PN) = T° ai ^ ; (14) 



Mo— 



J dM q 3 (p ,M) A(M,p N ) 



Mr, 



a/ A(/i, M, m^) y/X{p ,M, m„ 

q{fi, M) = , q (po, M) - 



2p ■ 2 p 

Here X(x, y, z) = (x 2 — (y — z) 2 )(x 2 — (y + z) 2 ) while is the pion mass; po, r° are the vacuum 
pole mass and width of the a\ spectral function with T° ai taken as 400 MeV. 

According to Eq. (fl4l) the a± width depends on the p meson spectral function A(M,pn) 
and, thus, is sensitive to the p meson in-medium effects. We stress, that such a 'coupling' of 
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FIG. 3: (Color online) The spectral function for the a\ meson for />jy = (solid line) and in case of 
the 'collisional broadening' scenario (dashed line) and the 'dropping mass + collisional broadening' 
scenario (short dashed line) for nuclear density pn = 3po. 

the di and p spectral functions has to be accounted for when drawing conclusions about an ex- 
perimental observation of chiral symmetry restoration by measuring the a± meson properties in 
heavy- ion experiments (e.g. the corresponding results in Ref. 40j should be strongly effected). 

Figure [3] shows the spectral function A ai (p, pn) for the a± meson for p^ = (solid line) 
and in the case of the 'collisional broadening' scenario (dashed line) and the 'dropping mass 
+ collisional broadening' scenario (short dashed line) for nuclear density p^ = 3po- Here the 
in-medium scenarios are related to both the a\ as well as the p meson. 

The production cross section a(s,M) for the tt + p —>■ a\ reaction (for all possible charge 
combinations) is calculated as follows, 

6tt 2 



(s,M,p N ) 



r ai 0, p N )A ai (s, p N ), 



(15) 



q 2 {s,M) 

where s = p? and the width T ai (s,piy) and momentum q(s, M) are defined by Eq. ([1] 

In Fig. H] we display the production cross section a(s, M) for the 7r + p — ► ai reaction 
for Piv = (solid lines) and in case of the 'collisional broadening' scenario (dashed lines) and 
the 'dropping mass + collisional broadening' scenario (short dashed lines) for nuclear density 
Pn = 3po- The upper part corresponds to the p meson mass of M = 0.77 GeV while the lower 
part shows the results for M = 0.3 GeV. As seen from Fig. H] the a\ production cross section 
in the 7r + p reaction is very sensitive to the in-medium scenarios as well as to the mass M of 
the initial p meson (even in the free case - cf. solid lines). 
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FIG. 4: (Color online) The production cross section a(s,M) for the n + p —* a\ reaction for p^ = 
(solid lines) and in the case of the 'collisional broadening' scenario (dashed lines) and the 'dropping 
mass + collisional broadening' scenario (short dashed lines) for nuclear density pn = 3po- The upper 
part (a) corresponds to M = 0.77 GeV while the lower part (b) shows the results for M = 0.3 GeV. 

D. Direct photon production from meson- meson bremsstrahlung 



The implementation of photon bremsstrahlung from hadronic reactions in transport ap- 
proaches is based on the 'soft photon' approximation. The soft-photon approximation (SPA) 



41] relies on the assumption that the radiation from internal lines is negligible and the strong 



interaction vertex is on-shell. In this case the strong interaction part and the electromagnetic 
part can be separated, so the soft-photon cross section for the reaction 1 + 2^1 + 2 + 7 can 
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be written as 

<i 3 cr 7 _ a a{s) 
d 3 q An 

s - ( Mi + M 2 f 
2Mf 

where Mi is the mass of the charged accelerated particle, M 2 is the mass of the second particle; 
qo,q are the energy and momentum of the photon. In (jl6p a(s) is the on-shell elastic cross 



?o— = —~g ( 16 ) 

_ s — \ivi\ + ivi 2 ) , . 

< s ) = a ( s )> 



section 
in Ref . 



or the reaction 1 + 2 — > 1 + 2. This approximation has also been employed by Haglin 



23| 



In Ref. 28] the photon bremsstrahlung from 7r + 7r and 7r + K elastic collisions has been 
calculated using the C/ em (l)-gauged meson-exchange model which includes the photon coupling 
to pseudoscalar and vector mesons. Indeed, such calculations go beyond the SPA model. How- 



ever, the direct comparison of SPA and the U em {l) models - cf. Fig. 4 in Ref. [28j| - show a 
very good agreement between the models, which allows us to use the simplified SPA formula 
for our purpose here. Thus, we have calculated the photon bremsstrahlung from all elastic 
meson-meson scattering mi + m<i — > mi + mi + 7 (where m = n, rj, K, K, K°, K*, K*, K*°), 
which occur during the heavy-ion collisions by applying the SPA formula (1161) . 

III. RESULTS FOR THE PHOTON PRODUCTION IN A+A COLLISIONS AT SPS 
ENERGIES 

Now we step on to the description of direct photon production from heavy-ion collisions ap- 
plying the HSD transport model incorporating all photon production channels from elementary 
reactions as described in the previous Section. 

In Fig. [5] we start with a comparison of the HSD results for the inclusive photon transverse 
momentum distribution for 10% central 158 A GeV 208 pb+ 208 Pb collisions with the experi- 
mental data from the WA98 Collaboration [15| (solid dots). We select photons in the pseudo- 



rapidity interval 2.35 < i] < 2.95 in the laboratory frame which corresponds to mid-rapidity 
in the center-of-mass frame. Our calculations show that such cut in pseudorapidity reduces 
the photon yield by about of factor 5, however, practically doesn't change the shape of the pt 
spectra. The WA98 data correspond to the transverse momentum interval 0.5 < px < 4 GeV/c 
while the HSD calculation includes also very low p? momenta but extends up to 2.5-3 GeV/c 
only since the high p? tail is very hardly reachable due to limitied statistics. 
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FIG. 5: (Color online) The inclusive photon transverse momentum distribution for central 158 A GeV 
208 Pb+ 208 Pb collisions at 2.35 < r] < 2.95. The solid symbols and arrows corresponds to the WA98 
data [if]], the lines correspond to the HSD results for the various channels as indicated in the legend 
(without medium effects for the vector mesons). 



Fig. presents the HSD calculations without medium effects for the vector mesons (the 
definition of the lines are indicated in the legend). As seen from Fig. [5] the dominant channels 
for the inclusive photon production are the photon from 7r and rj decays (tt — > 7 + 7 and 
7] — > 7 + 7), whereas other channels are down by more than order of magnitude. The HSD 
results agree very well with the experimental data which is, indeed, expectable since the HSD 
model provides a good description of the pion transverse momentum spectra at SPS energies 



42j ] and predicts a meson scaling [43]. 

In order to obtain the information about the direct photon production, one needs to subtract 
the 'background' contributions which are dominated by the mesonic decay processes. While the 
7T° and rj spectra are measured directly by the same experiments [3], it is possible to estimate 
their photon decay in a reliable way. Indeed, the life times of 7r's and r/'s are large such that 
they basically decay at the end of reaction stage, i.e. in vacuum. However, the situation with 
the rf — > p + 7, <jj — > 7T° + 7, cf) — > rj + 7, d\ — >• 7r + 7 decays are not so transparent since these 
mesons cannot be mesured directly by the WA98 detector. 
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FIG. 6: (Color online) The nucleon density distribution at the photon emission point (in units of 
p ) for central 158 A GeV 208 Pb+ 208 Pb collisions at 2.35 < rj < 2.95. The definition of the lines is 
indicated in the legend. 

The WA98 Collaboration has subtracted the contribution of hadronic decays by calculations 
based on the assumption of m^-scaling with the same slope of the m^-spectra as measured 
in the 7r° spectrum and with relative normalizations Rh a dron/n° (equivalent to the asymptotic 
ratios for px — > oo) [l5]. However, such a procedure is a rought approximation especially for 
the short living resonances which can emit photons inside the hot fireball shortly after creation 
and, thus, can not be reconstruct in experiment. This statement is illustrated in Fig. [H] which 
shows the density distribution at the photon emission points from rj', u, <fr, a>\ decays for central 
158 A GeV 208 Pb+ 208 Pb collisions at 2.35 < rj < 2.95. One can see from Fig. 1 that there 
is some fraction of LU,(p,ai mesons that emit photons for > Po/2 (e.g. ~ 25% of cj's and 
~ 70% of cti's), whereas a relatively long leaving rj' meson decays dominantly at low densities. 

In order to demonstrate the possible contribution of the in-medium decays of t]', u, <fr, a\ 
mesons to the direct photon spectra reported by the WA98 Collaboration we show in the r.h.s 
of Fig. [7] the invariant mid-rapidity spectra including the contribution of the rf, u, <f), a± photon 
decay only for p^ > p /2 whereas the l.h.s. of Fig. [7] shows the total contribution of all 
hadronic channels at all densities (for comparison). Indeed, such a cut in density leads to 
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FIG. 7: (Color online) The invariant direct photon multiplicity for central 158 A GeV 208 pb+ 208 Pb 
collisions at 2.35 < rj < 2.95. The solid symbols and arrows corresponds to the WA98 data fl5|. 
the lines stand for the HSD calculations without (upper row - (a),(c)) and with (lower row - (b),(d)) 
medium effects (i.e. collisional broadening) for the vector mesons. The definition of the lines is 
indicated in the legend. The solid line for px > 1-5 GeV/c gives the contribution from the prompt 
photons. The right panel (c),(d) shows the spectra including the contribution of the r( , u>, cj>, a\ photon 
decay only for > po/2 whereas the left panel (a),(b) shows the total contribution of all hadronic 
channels at all densities. 



a substantial reduction of the hadronic decay 'background' especially for low px, however, it 
becomes essential at pt < 0.7 GeV/c. Thus, our calculations provide the scale of possible 
uncertainties in the experimental background subtraction which is important for the physical 
interpretation of the experimental results. 

One can see from Fig. [7] that the dominant channel for low pt phot on p roduction is meson- 
meson bremsstrahlung. This is in line with the analysis by Haglin 27J] and Liu and Rapp 
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FIG. 8: (Color online) The invariant direct photon multiplicity for central 158 A GeV 2U8 pb+ 2U8 Pb 



collisions at 2.35 < r\ < 2.95. The solid symbols and arrows corresponds to the WA98 data 
the lines stand for the HSD calculations: the upper solid line with the stars shows the sum over all 
channels, the solid line corresponds to the sum over all meson-meson bremsstrahlung channels which 
are specified by the various lines. The definition of the lines is indicated in the legend. 



281 ] . The contributions of the other channels such as direct photon production by meson- 



meson collisions n + n — > p + 7 and tc + p — > n + 7 are found to be suppressed relative to 
the meson-meson bremsstrahlung. Also we have investigated the influence of in-medium effects 
such as collisional broadening for the vector meson spectral functions which leads (e.g.) to 
enhancement of low mass dilepton production 17|, ll8|]. However, we here obtain only a small 
impact on the direct photon production which is hidden below the dominant bremsstrahlung 
contribution. The same holds also for other in-medium scenario such as a dropping mass of the 
vector meson. 

In Fig. [8] we specify the meson-meson bremsstrahlung by showing the contributions of 
different sub-channels 7r + 7r^7r + 7r-|-7, ir + r]—>-7r + r] + '-f, tc + K^tt + K + 'j with 
K = (K + , K~), 7r + K° -> 7T + K° + 7, vr + K* -> 7T + K* + 7 and vr + K*° -> vr + K*° + 7. 
We stress again that only the elastic meson-meson collisions have been accounted here. Indeed, 
the final result is very sensitive to the meson-meson elastic cross section which is, however, 
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very hard to measure experimentally. We have used a 10 mb elastic cross section for all meson- 
meson scatterings in the present study. As follows from Fig. [8] the dominant contribution 
stems from tt + 7r elastic scattering; the contributions from ir + K*°, ir + K° and tt + rj are also 



visible at low px- This is again consistent with Refs. [27|, |28( where the authors found that the 
photon emission rate from pion-kaon scattering is at the ~ 20% level of the one from pion-pion 
scattering. 

We mention that the solid line at high px in Fig. [7] shows the estimated contribution from 



the prompt photons. For that we have convoluted the pp data fit from Ref. |44j with the nuclear 
overlap function Tab- Indeed, it gives a lower estimate for the prompt photons since the nuclear 
effects (as e.g. Cronin effect) have been ignored. For futher discussions on the contributions of 



the prompt as well as thermal photons we address to Refs. 26|, |44j and references therein 



IV. RESULTS FOR THE PHOTON PRODUCTION IN P+A COLLISIONS AT SPS 
ENERGIES 

In order to draw solid conclusions on the direct photons from heavy-ion collisions one needs 
some 'reference systems', i.e. to compare to p+p and/or p+A data. Recently the WA98 
Collaboration has provided preliminary data on p+C and p+Pb collisions at ^/s = 17 .4 GeV. 
In this Section we present the HSD results for these systems. 

We start again with the inclusive invariant photon spectra from p+C and p+Pb collisions 
at mid-rapidity for 160 A GeV which are depicted in Fig. [9j The solid dots corresponds to the 
preliminary WA98 data [3], the various lines stand for the HSD calculations. Again, as in a 
case of Pb+Pb, the dominant contribution comes from the tt and r\ photonic decays. As seen 
from Fig. [9] HSD provides a very good description of the inclusive experimental data. 

In order to focus on the direct photon contribution the WA98 collaboration has subtracted 
the contributions of the hadronic decays in a similar way as for the Pb+Pb collisions. However, 
due to experimental uncertainties in the 'background ' subtraction the preliminary WA98 data 
correspond to the upper limits on direct photon production [19( and are shown by the (down) 
arrows in Fig. [10] for p+C (upper part) and p+Pb (lower part). The various lines in Fig. [10] 
stand for the HSD result. Here we do not apply any cuts in density in order to obtain an upper 
estimate for the hadronic decays. Contrary to the Pb+Pb collisions the contribution from the 
meson-meson bremsstrahlung is suppressed compared to the rj' and uj decay due to the small 
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FIG. 9: (Color online) The inclusive invariant photon multiplicity for 158 A GeV p+ 12 C (upper part 
- (a)) and p+ 208 Pb (lower part - (b)) at 2.3 < r\ < 3.0. The solid dots corresponds to the preliminary 
WA98 data [19[; the lines stand for the HSD calculations. The definition of the lines is indicated in 
the legend. 

meson density in p+A collisions. As seen from Fig. \TU\ the HSD results are much below the 
upper limits from WA98. Since it is unlikely that there are some new channels for photon 
production in p+A reactions, which are not visible in Pb+Pb collisions, we expect that the 
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FIG. 10: (Color online) The invariant direct photon multiplicity for 158 A GeV p+ 12 C (upper part - 
(a)) and p+ 208 Pb (lower part - (b)) at 2.3 < r/ < 3.0. The solid dots corresponds to the preliminary 



WA98 data 
the legend. 



19( |; the lines stand for the HSD calculations. The definition of the lines is indicated in 



final data might go down, too. Thus, it is quite important to obtain precise p+A data from 
the experimental side in order to check the theoretical models more accurately. 
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V. SUMMARY 



In this study a detailed analysis of the low px photon production in p+C, p+Pb and Pb+Pb 
collisions at 160 A GeV has been presented within the microscopic HSD transport approach 
that incorporates a full off-shell propagation of the vector mesons [l7|. We stress, that the HSD 
approach contains only hadronic degrees of freedom and strings and doesn't include the phase 
transition from a QGP to hadronic matter explicitly; the partonic interactions are treated in 
HSD only inside the strings (and in terms of leading quarks and diquarks). However, the model 
has all non-equilibrium dynamical features and demonstrates a good ability in describing a 
variety observables for heavy-ion, proton- and pion- induced reactions. Accordingly, we apply 
our model to study the direct photon production from hadronic sources which dominate at low 
p T . 

In particular the following hadronic sources for direct photon production in heavy-ion colli- 
sions have been incorporated: 

- photon emission from elementary meson-meson rescatterings, where the processes irp — > 7r7 
and 7T7T — > p7 are dominant due to the high pion production rate. The novel issue here is that 

p meson and extended the vacuum cross 



29| for the in-medium case with the full 



we accounted for the off-shellness of the initial/fina 
sections for these processes from Kapusta et al. 
off-shell p meson spectral function. It allows us to investigate the influence of the in-medium 
effects such as 'collisional broadening' on the photon pt spectrum. 

- meson-meson bremsstrahlung from the elastic meson-meson scattering m\ + m 2 — > mi + 
m 2 + 7 (where m = ir, rj, K, K, K°, K*, K*, K*°), which we accounted for by applying the SPA 
formula ([IS]) . 

We have found that the enhancement of the low pt photon emission from elementary meson- 
meson rescattering due to a collisional broadening of the vector-meson spectral functions is 
hardly visible in the final spectra which are dominated by bremsstrahlung type processes. 



27 



28] 



Thus our non-equilibrium dynamical calculations support the early findings in Refs. 
based on the hadron gas models. 

Also we have investigated the uncertainties in the extraction of the direct photon yield 
from the measured inclusive photon spectra which are dominated by the hadronic decays, 
in particular 7r°, rj, rj', u and a\. Here the HSD model shows a very good agreement with 
the measured inclusive photon transverse momentum distribution. While the 'background' 
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contributions from the 7r° and rj decay can be subtracted in a reliable way - since these mesons 
are measured by the same WA98 experiment - there remains a problem in the subtraction of 
the contributions from short living resonances which decay inside the hot and dense fireball. 
We have estimated the possible contribution of such in-medium decay processes by comparing 
two calculations - in the first we have selected the photons coming from the mesonic decays 
inside the fireball by applying a density cut pn > Po/2, in the second we have accounted for 
the mesonic decays over the full time history. We found that the density cut reduces drastically 
the final yield which depends directly on the actual value of the cut. However, some amount 
of photons from the inside oo decay for p^ > po/2 is visible in the final spectra and might be 
misidentified as direct photons. 

For reference we have also studied the photon production for p+C and p+Pb collisions at 
160 GeV. Again the HSD model provides a good description of the preliminary experimental 
data from the WA98 Collaboration on the inclusive spectra, however, is much below the upper 
limits for the direct photons (even without a subtraction of r)',u,(j> and a x photon decays). In 
this respect it will be very helpful to have the final experimental data (with high accuracy). 
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